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ABSTRACT 

With  high-speed  digital  signal  processing  (DSP)  and  signal  digitization  hardware  becoming 
available  at  a reasonable  cost,  it  is  now  practical  to  propose  DSP  demodulation  for  future  designs. 
This  paper  presents  the  necessary  theory  and  governing  equations  for  practical  DSP  demodulators 
that  can  be  implemented  for  a wide  variety  of  applications.  The  Generalized  Costas  Demodulator 
(GCD)  will  be  emphasized  because  it  is  capable  of  recovering  the  demodulation  functions  of  an 
arbitrary  bandlimited  signal. 


1.  INTRODUCTION 

The  use  of  digital  signal  processing  for  modem  designs  has  been  very  popular  the  past  few  years. 
The  effective  application  of  DSP  to  these  designs  requires  a good  understanding  of  the  governing 
equations  associated  with  modulation/demodulation  and  the  signal  digitization  process.  In  particular, 
a good  understanding  of  quadrature  detectors,  aliasing,  decimation,  and  the  process  of  carrier 
recovery  is  required.  This  paper  has  been  written  to  contribute  to  this  understanding  by  presenting 
the  necessary  governing  equations. 

A good  DSP  design  depends  on  a good  signal  digitization  technique.  The  designer  must  be 
concerned  with  the  linearity,  dynamic  range,  and  noise  performance  of  the  digitizer.  It  is  especially 
difficult  to  achieve  the  dynamic  range  needed  for  most  receiver  designs  because  of  1 ) digitizer  range 
limitations  and  2)  quantization  noise  produced  during  digitization.  Usually,  automatic  gain  control 
(AGC)  or  hard  limiters  are  used  to  reduce  the  dynamic  range  requirements  for  the  digitizer.  Recently 
[1]  [2]  sigma-delta  modulation  has  been  proposed  for  use  in  signal  digitizers  as  a method  of 
achieving  high  dynamic  range  and  good  linearity. 

The  governing  equations  for  modulation/demodulation  will  now  be  presented.  It  will  be  assumed 
that  adequate  signal  digitization  has  been  achieved.  An  illustrative  modulation/demodulation 
example  is  included  in  the  appendix. 


2.  CONCEPTUAL  SIGNAL  MODEL  - TRANSMITTER 


To  adequately  describe  DSP  demodulation,  it  is  first  necessary  to  develop  precise  conceptual 
models  for  the  transmitted  and  received  signals.  The  fundamental  model  underpinning  all  analysis 
is  that  of  a noiseless  narrowband  signal  that  can  be  analytically  represented  in: 

1 . amplitude/phase  form, 

dt)  = di)  Accos[2nfct+  <}>c+  <()(d]  (l) 

2.  or  quadrature  form  [3], 

dt)  = a(t)Accos,[2nfct+  <|>c]  - bit)  A c sin[27t  fct+  ((>  f) 

3.  or  upper/lower  sideband  form  [4]: 

sit)  = [uii)  + 1{£)  + A ccob[2k  fct+ § ^ - [uit)  - ~lii)]  A csm[2n  fct+ § ^ 

The  signal  carrier  is  modeled  as  a cosine  with  amplitude  Ac,  frequency  fc,  and  constant  phase,  4>c- 
The  constant  phase  will  be  used  to  model  the  effect  of  an  arbitrary  reference  time  at  the  receiver. 
The  effect  of  modulation  on  the  carrier  is  modeled  using  the  envelope  and  phase  functions,  i(t)  and 
<t>  it) , the  quadrature  functions,  ait)  and  bit) , or  the  upper  and  lower  sideband  functions,  u(t)  and  At) . 
These  modulation  functions  all  represent  some  coded  form  of  an  analog  or  digital  information 
source.  For  this  narrowband  model,  the  modulation  of  the  carrier  must  be  restricted  such  that  the 
spectrum  of  dt)  is  zero  outside  the  frequency  range  fY  to  f2  where  if^-f)  « fc  and  all  the 
transmitted  spectrum  is  passed  by  the  preselector  and  IF  filters  in  the  receiver.  An  unmodulated 
carrier  will  be  represented  by  the  constants  dt)  = 1 and  4>(d  = 0. 


Using  trigonometric  identities  [5],  the  modulation  functions  can  be  related  as  follows: 
ait)  = r<£)cos[4)(d]  bit)  = Xdsin[(|)(£)] 

bit) 


At)  = \Ja2it)  + b2it)  4>(d  = arctan 


ait) 
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The  modulation  functions  may  assume  any  real  value  except  for  the  envelope  function,  z{i) , which 
must  always  be  non-negative. 

The  amplitude  spectrum  of  dt) , if  it  exists,  can  be  obtained  by  applying  the  Fourier  transform 
to  the  quadrature  form  of  the  signal.  Because  sit)  is  a narrowband  signal,  the  spectrum  is  nonzero 
only  near  the  frequencies  ± /.  Furthermore,  a sinusoidal  carrier  frequency  component  is  present  only 

if  A{  /=/)  or  B(  it/)  contains  a Dirac  delta  function  [3].  The  tilde  (~)  denotes  the  Hilbert  transform 
needed  to  generate  SSB  signals.  The  constant,  kc,  represents  the  magnitude  of  any  carrier  that  is 
transmitted  with  the  two  sidebands. 

3.  CONCEPTUAL  SIGNAL  MODEL  - RECEIVER 

The  simplified  signal  model  at  the  receiver  will  include  the  effects  of  additive  channel  noise  and 
carrier  phase  shift  due  to  propagation  delay.  Doppler,  dispersion,  slow  carrier  amplitude  variations, 
and  multipath  effects  will  not  be  modeled  in  this  paper.  The  constant  phase.  4>c,  represents  the 

carrier  phase  shift  observed  at  the  receiver  due  to  propagation  delay  and  an  arbitrary,  but  known, 
local  reference  time. 

Additive  channel  noise,  nit) , can  be  thought  of  as  an  additional  source  of  modulation  that  is 
random  and  unwanted.  The  composite  signal,  sc(t) , in  a narrowband  IF  filter  is  the  sum  of  the 
noiseless  signal  and  narrowband  noise: 

sc(t)  = i(t)  Accos[2%fct+  4>c+  4>  (/)]  + nit)  (8) 

Since  the  channel  noise  is  narrowband,  it  can  be  represented  in  quadrature  form  [6]  as: 

lit)  = nc(t)co&[2^fct+  (f)c]  - 73s(z)sin[2ji/r+  (9) 

The  quadrature  noise  components,  nc(t)  and  ns(t) , are  baseband,  zero-mean,  uncorrelated,  Gaussian 
noise  signals  with  a spectral  shape  that  is  the  same  as  the  IF  filter. 
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The  quadrature  form  of  the  composite  IF  signal  is  easily  computed  as  the  sum  of  the  quadrature 
forms  of  the  modulated  carrier  and  the  narrowband  noise, 


sc(t)  = a(t)  Accos[2n  fct+  (J>cl  - bit)  Acsin[2n  fct+  (J>cl 

(10) 

where: 

ait)  = a(i)  + nc(i)IAc  = X/)cos[(j)(d] 

(11) 

bit)  = bit)  + nJit)/Ac  = Xz)sin[(j)(z)] 

(12) 

and  the  "hat"  indicates  that  the  function  has  been  corrupted  by  noise.  The  relationships  among  all 
the  noise-corrupted  modulation  functions  are  the  same  as  noiseless  ones  given  in  (4)  through  (7). 


It  is  important  to  note  that  the  channel  noise  terms  are  always  scaled  by  the  amplitude  of  the 
carrier  signal,  i.e.  the  representation  chosen  for  this  model  accounts  for  the  signal-to-noise  ratio  in 
a "natural"  way.  The  computation  of  the  sideband  signals  requires  the  implementation  of  a Hilbert 
transform  in  the  receiver. 


4.  THE  DEMODULATION  PROCESS 

Demodulation  is  the  process  of  recovering  the  noise-corrupted  modulation  functions  from  the 
received  signal.  It  can  also  be  thought  of  as  carrier  removal  or  "carrier  stripping".  To  properly 
demodulate  a signal,  the  carrier  parameters  (amplitude,  frequency,  and  phase  constant)  must  be 

known  or  estimated  by  the  receiver.  In  other  words,  {aii),tii)\,  {)(d,  £riT)},  or  {i(d,<i>(d}  cannot,  in 
general,  be  extracted  from  the  carrier  without  a knowledge  of  the  carrier  parameters.  The  carrier 
amplitude  governs  the  linearity  and  dynamic  range  requirements  of  the  receiver  as  well  as  the  design 
of  the  signal  digitizer.  Carrier  frequency  changes  (not  those  due  to  frequency  modulation)  and/or 
frequency  channelization  require  that  the  receiver  be  precisely  tuned.  Coherent 
modulation/demodulation  schemes  also  require  precise  estimation  (tracking)  of  both  carrier 
frequency  and  carrier  phase. 

The  noise-corrupted  modulation  functions  will  be  called  demodulation  functions  and  Ac,  fc  and  4>  c 

will  be  called  carrier  parameters.  To  be  effective,  a receiver  must  produce  good  estimates  of  the 
carrier  parameters  and  then  successfully  extract  the  demodulation  functions.  However,  the  original 
modulation  functions  can  never  be  exactly  recovered  because  they  have  been  irretrievably  corrupted 
by  noise. 
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The  most  general  approach  to  demodulation  is  to  implement  "I-Q"  or  quadrature  detection  and 
then  remove  the  carrier  by  obtaining  good  estimates  of  the  carrier  parameters.  Figure  1 shows  a 
basic  quadrature  detector.  To  begin,  the  received  signal  is  multiplied  by  the  local-carrier  quadrature 
functions  as  shown  in  the  following  equations: 

x^t)  = se(d2^4ocos[0o(4)]  (13) 

y^t)  = sc(  2)  2 ^4Q  sin  [0  0(  *)  ] (14) 

where: 

0o(d  = 2nf0t+  (J)0-  r)(d  (15) 

—Q0(t)  = 2nf0  - —r\(t)  (16) 

dt  0 dt 

The  parameters  Aq,  fQ,  <t>0,  and  0 0 (t)  for  the  local  carrier,  represent  estimates  of  Ac,  fc,  <t>  c,  and  Q c(i), 

respectively.  The  additive  noise  term,  r)  (t) , is  the  error  in  the  receiver's  estimation  of  the  carrier 
phase  function. 

When  the  trigonometric  products  representing  the  mixer  outputs  are  expanded  using  trig 
identities,  the  equations  representing  the  mixer  outputs  can  be  expressed  as: 

x^t)  = i(£)A'cos[0c(£)  + 0o(d  + <i>(2)]  + X£)iTocos[0c(£)  - QQ(t)  + $(0]  (17) 

and, 

yx(t)  = i-(d.K'osin[0c(f)  + 0Q(d  + 4>(dl  - Xd  Wosin[0c(z)  - 0o(f)  + 4>(z)]  (18) 

Lowpass  zonal  filtering  [7]  (an  "ideal"  lowpass  filter)  can  be  applied  to  the  multiplier  outputs  to 
reject  the  "sum"  terms  and  pass  the  "difference"  terms  so  that  the  filter  outputs  are, 

x2(t)  = r(t) K0cos[Q c(t)  - 0o(d  + $(£)]  = r(t)KQ costy  (2)]  (19) 

y2(t)  = - Xt)  Wosin[0c(d  - 0o(d  + <f>(dl  = - XdA"0 sinter  (dl  (20) 

where  the  substitution,  (t)  = Qc(t)  - Q0(t)  + <f>  (2),  will  be  used  to  simplify  the  notation. 


This  widely-known  result  is  valid  only  if  the  received  signal  is  narrowband  and  the  zonal  filter 
passes  an  undistorted  version  of  the  entire  spectrum  of  the  difference  term.  In  practical  designs,  the 
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zonal  filter  is  implemented  using  a lowpass  filter  having  a stopband  attenuation  that  adequately 
suppresses  the  sum-frequency  output  while  still  passing  the  entire  frequency  spectrum  of  the 
difference-frequency  output.  If  the  lowpass  filter  distorts  the  difference-frequency  spectrum,  this 
result  assuming  zonal  filtering  is  not  valid.  The  technical  literature  frequently  uses  the  symbols  I and 
Q to  represent  these  quadrature  outputs.  The  quadrature  outputs  of  the  lowpass  zonal  filters  contain 
all  of  the  recoverable  information  about  the  transmitted  signal  as  well  as  the  undesired  noise  caused 
by  both  the  channel  and  the  imperfect  estimation  of  the  local-carrier  phase  function. 

The  quadrature  detector  outputs  can  be  processed  using  summing,  scaling,  squaring,  square  root, 
inverse  tangent,  and  Hilbert  transform  operations  to  estimate  all  of  the  modulation  functions  as 
follows: 


rit)  — rit)  = — \jx2it)  + y2it) 


(21) 


c()  it)  — cf>  it)  — 4*  (d  = arctan 


~T2(4> 

x2it) 


(22) 


ait)  - ait)  — 


bit)  - bit)  — 


lit)  - lit)  - hx2it)-y2it)] 

2 


(23) 


(24) 


(25) 


uit)  - uit)  - — [x2ii)  + y2ii>]  (26) 

2 

The  estimates  for  the  phase  function,  quadrature  functions,  and  the  upper  and  lower  sideband  signals 
are  valid  only  if  the  receiver  is  "tracking"  the  carrier,  i.e.  only  if  0O(£)  + q it)  = 0c(d-  The  double 

approximations  are  used  to  indicate  the  separate  influences  of  channel  noise  and  local  oscillator 
noise  in  estimating  all  the  outputs  that  are  dependent  on  phase. 

For  specific  types  of  modulation,  there  can  be  ways  to  recover  the  demodulation  functions  that 
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avoid  the  need  for  a quadrature  detector.  For  example,  a square-law-based  envelope  detector, 
represented  by  the  following  equation,  can  be  used  to  recover  the  envelope  function  without  the  use 
of  a quadrature  detector: 


= 


hXt ) * 


sHt) 


-Arr(t) 
2 c 


For  phase  modulation,  a simple  phase-locked  loop  can  be  used.  These  two  methods  are  represented 
by  the  block  diagram  of  Figure  2.  In  this  figure,  x5  is  the  output  of  the  bandpass  limiter  which 
removes  the  envelope  modulation  and  preserves  the  phase  modulation.  The  PLL  phase  output  is  x6 
which  is  an  estimate  of  the  phase  modulation  function.  The  lower  branch  of  the  figure  represents 
the  processing  needed  for  envelope  demodulation. 


In  summary,  a generalized  demodulator  must:  1)  have  a prior  knowledge  of,  or  estimate,  the 
carrier  parameters,  2)  generate  the  estimated  carrier  for  use  in  an  I/Q  detector  (or  a correlator),  and 
3)  recover  the  demodulation  functions. 


(27) 


5.  GENERALIZED  COSTAS  DEMODULATOR 

The  Generalized  Costas  Demodulator  (GCD)  of  Figure  3 represents  a general  method  of  signal 
demodulation.  Its  essential  elements  are,  1)  a quadrature  detector,  2)  carrier  estimation,  and  3)  a 
Hilbert  transform.  The  signal  processing  blocks  to  the  right  of  the  quadrature  detector  implement 
Equations  21  through  26. 

This  general-purpose  demodulator  block  diagram  represents  all  of  the  governing  equations 
necessary  to  demodulate  any  type  of  signal.  It  can  be  implemented  using  DSP  if  proper  bandpass 
signal  digitization  is  used.  The  next  section  will  discuss  various  issues  for  this  digital 
implementation. 


6.  DISCRETE-TIME  DEMODULATION 

Discrete-Time  (or  DSP)  demodulation  refers  the  concept  of  implementing  the  demodulation 
process  in  digital  signal  processing  hardware  and  software  [8]  [9]  [10]  [1 1]  [12].  A successful  DSP 
demodulator  design  requires  the  proper  use  of  1)  sampler  prefiltering,  2)  sampling  design  (including 
aliased  sampling),  3)  digital  lowpass  filtering,  and  4)  decimation.  In  addition,  adequate  algorithms 
must  be  available  for  the  following  operations:  1)  multiplication,  2)  squaring,  3)  square  root,  and  4) 
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arctangent.  Care  must  be  taken  to  insure  that  the  noise  caused  by  number  truncation  in  the 
algorithms  does  not  degrade  the  signal-to-noise  ratio  in  the  demodulator. 

The  signal  can  be  sampled  1)  after  the  bandpass  filter,  2)  after  the  mixers,  or  3)  after  the  lowpass 
filters.  A knowledge  of  the  signal  spectrum  at  each  point  in  the  block  diagram  is  necessary  if  the 
design  is  to  avoid  signal  aliasing.  Also,  this  knowledge  will  help  the  designer  to  implement 
decimation  where  it  will  be  the  most  effective. 

Sampling  after  the  Bandpass  Filter.  Baseband  or  aliased  passband  sampling  [10]  can  be  used  after 
the  bandpass  filter  as  long  as  the  sampling  frequency,  fs,  is  chosen  such  that  the  spectrum  of  the 
"sum"  frequency  term  can  be  adequately  filtered  while  still  passing  the  "difference"  frequency  term. 
This  requirement  is  imposed  because  of  the  multiplication  that  is  performed  in  the  sine  and  cosine 
mixers.  When  using  the  governing  equations,  the  carrier  frequency,  fc,  must  be  interpreted  as  the 
carrier  frequency  after  aliased  sampling.  In  actual  designs,  some  spectral  folding  is  acceptable  as 
long  as  the  lowpass  filter  can  still  effectively  remove  all  of  the  spectrum  associated  with  the  "sum" 
frequency,  i.e.,  it  performs  like  a zonal  filter.  Decimation  can  be  applied  to  the  outputs  of  the 
lowpass  filters  to  reduce  the  processing  requirements  for  subsequent  DSP  operations.  Care  must  be 
taken  not  to  decimate  to  the  point  where  subsequent  operations  might  produce  undesirable  spectral 
folding. 

Sampling  after  the  Mixers.  Baseband  sampling  can  be  used  after  the  mixers  as  long  as  the  folding 
frequency,  fs/2,  is  high  enough  to  prevent  the  folding  of  the  "sum"  frequency  terms  into  the 
equivalent  passband  of  the  lowpass  filters.  Some  decimation  can  be  applied  to  the  outputs  of  the 
lowpass  filters  but  the  number  rate  should  still  be  at  least  four  times  the  cutoff  frequency. 

Sampling  after  the  Lowpass  Filters.  If  the  carrier  frequency  is  too  high  for  adequate-quality  aliased 
sampling,  it  is  still  possible  to  sample  after  the  lowpass  filters  and  then  use  DSP  algorithms  for  the 
remaining  demodulation  processing.  Because  nonlinear  operations  such  a squaring  are  used,  the 
sampling  rate  (number  rate)  at  the  filter  output  must  be  much  higher  that  the  Nyquist  rate.  For 
example,  Fourier  analysis  shows  that  the  spectrum  of  x2  or  y2  will  be  doubled  after  squaring  so  the 
number  rate  must  be  at  least  four  times  the  cutoff  frequency  of  the  lowpass  filter.  Designers  should 
carefully  analyze  the  output  spectrum  for  each  specific  design  and  choose  the  number  rate 
accordingly. 

To  summarize,  the  signal  spectrum  at  the  sampling  point  and  all  subsequent  numerical  operations 
should  be  carefully  studied  before  selecting  the  sampling  rate  and  the  amount  of  decimation  to  be 
applied.  Many  designs  can  be  effective  with  aliased  sampling  of  the  bandpass  filter  output  and  a 
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number  rate  at  the  quadrature  outputs  of  at  least  four  times  the  lowpass  filter  cutoff  frequency. 


7.  CARRIER  ESTIMATION 

The  generation  of  a local  replica  for  the  carrier  is  required  for  proper  operation  of  the 
Generalized  Costas  Demodulator  or  any  synchronous  demodulator.  In  order  to  generate  a local 
replica,  the  receiver  must  be  capable  of  estimating  the  carrier  parameters  Ac,  /.  and  4>  c by  using  one 

or  more  of  the  signals  represented  by  equations  10,  19,  or  20.  The  required  accuracy  of  the 
estimation  for  each  parameter  will  depend  upon  the  modulation  type  and  the  available  signal-to-noise 
ratio. 

Carrier  estimation  is  accomplished  in  a variety  of  ways  using  a combination  of  filtering,  phase 
tracking,  and  nonlinear  operations  such  as  squaring  or  bandpass  limiting.  The  governing  equations 
are  summarized  below. 

Narrowband  filtering.  If  a pilot  carrier  or  other  carrier  frequency  component  is  present  and  the 
modulation  sidebands  are  sufficiently  separated  from  the  carrier,  narrowband  filtering  can  be  used: 

sc(t)*hn(i)  ~ Acco&[2n  fct+  <}>c]  (28) 

In  this  case,  the  carrier  phase  constant  also  models  the  delay  in  the  filter. 

Bandpass  limiter.  When  only  envelope  modulation  is  present,  a bandpass  limiter  (hard  limiter 
followed  by  a narrowband  filter)  can  be  used: 

sgn{sc(i)}  * hn(i)  = —cos[2nfct+  (J)c1  (29) 

The  signum  function,  sgn{ },  [3]  represents  the  action  of  the  hard  limiter.  The  symmetrical  limiter 
clipping  levels  are  ±a. 

Phase-locked  Loop.  For  some  types  of  modulation,  phase-locked  loop  can  be  used  to  regenerate  the 
carrier.  For  best  performance,  a PLL  is  often  used  in  conjunction  with  hard  limiters  and/or  filters. 

The  PLL  local  oscillator  produces  a waveform  with  a fundamental  frequency  component  of: 

s0(t ) = A0sin[2rt  fQt+  (|>0-  r\  (dl  (30) 

where  f0  = fc  and  <j>0=  <j>c  when  the  loop  is  properly  tracking. 
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Costas  Loop.  The  Costas  loop  is  similar  to  the  PLL  in  that  it  controls  a local  oscillator  that  tracks 
the  carrier  frequency  and  phase.  The  error  signal  for  the  loop  is  obtained  by  averaging  the  product 
of  the  quadrature  outputs  as  represented  below. 

(x2 {t) y2 (z)^  = - 2*2(z)  cos[i|f  (t) ] sin[i|f  ( t ) ] 

= — ^-(^2(^)sin[2i|f  (dl)  (31) 

~ --/£,(2i|f  (d)  ~ 0 

The  control  loop  works  to  keep  the  average  value  of  \Jj  (t)  at  zero. 


Squaring  Loop.  The  squaring  loop  [13]  is  used  for  binary  phase-shift  keyed  (BPSK)  signals  and 
works  by  filtering  or  phase  locking  to  the  second  harmonic  of  the  squared  narrowband  signal 
represented  below.  Notice  that  by  doubling  the  phase  function,  the  phase  modulation  has  been 
removed. 


cos[20c(d  + 2<j>(z)] 


(32) 


The  carrier  is  regenerated  using  a divide -by-two  on  the  PLL  oscillator  output. 


Multiplier,  m-th  order.  The  squaring  loop  idea  can  be  applied  to  m-ary  phase  modulation  by 
obtaining  the  m-th  harmonic  of  the  narrowband  signal  as  represented  by  the  equation: 


s™(t)  = rm(t)  cos^tB  (t)  + $ (2)] 


= r 


n(t) 


L COS  [ 2220  C)  + 222(f)  (z)] 


(33) 


Derivative.  An  estimate  of  the  carrier  frequency  can  be  obtained  by  applying  the  derivative  to  the 
output  of  the  bandpass  limiter  [14],  This  result  is: 


a. [ 
dtl 


2ti/+  — <f>  (t) 
dt 


[sg22j;Sc(z)}  * hn(t)^  = - — sin[27i  fct+  (f)c+  (2)] 

sin[27i  f t+  (f)  + (2)] 


8 af  + — — ^ (t) 
n dt 


(34) 


The  average  value  of  the  envelope  function  of  the  result  is  used  for  the  estimate: 


fc 


8 a 


— (-[8 


7t  dt 


(35) 


For  this  approach  to  be  valid,  the  average  value  of  the  derivative  of  the  phase  modulation  function 
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must  be  zero. 


For  some  types  of  modulation  it  is  not  possible  to  obtain  an  estimate  of  the  carrier  parameters 
directly  from  the  signal.  Single-sideband  modulation  is  an  example  of  this.  For  these  cases,  the 
receiver  must  have  an  accurate  prior  knowledge  of  the  carrier  frequency. 


8.  CONCLUSIONS 

The  successful  development  of  hardware  and  software  designs  for  DSP  demodulation  requires 
a good  understanding  of  the  governing  equations  for  demodulation  and  the  aliasing  problems 
associated  with  discrete-time  nonlinear  operations.  Using  the  governing  equations,  it  is  possible  to 
develop  one  design,  the  Generalized  Costas  Demodulator,  that  will  work  for  any  type  of  modulation. 
For  new  designs,  aliased  sampling  at  the  output  of  the  bandpass  prefilter  should  be  considered  first. 
A number  rate  at  least  twice  the  Nyquist  rate  should  be  used  at  the  output  of  the  lowpass  filters  in 
the  quadrature  detector. 


9.  APPENDIX 

An  example  will  now  be  given  to  illustrate  how  the  modulation  functions  are  selected  for  a well- 
known  modulation  type,  Minimum-Shift  Keying  (MSK).  In  MSK,  one  of  four  possible  waveforms 
are  sent  using  a phase-coherent  frequency-shift  keying  technique.  The  four  possible  waveforms  are 
generated  by  selecting  one  of  two  possible  frequencies  and  one  of  two  possible  phases.  The 
selections  are  restricted  such  that  either  the  slope  (MSK-1)  or  the  magnitude  (MSK-2)  of  the 
transmitted  signal  is  zero  at  the  time  the  next  waveform  is  selected. 

To  insure  these  smooth  transitions,  the  carrier  frequencies  are  restricted  to  the  values, 

1 

n+  — 

f _ 2 (36) 

where  n is  an  integer  (usually  large)  and  Tb  is  the  interval  between  data  bits.  The  two  frequencies 
transmitted  by  MSK  are  related  to  Tb  and  the  carrier  frequency  as  follows: 
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£ = -B-  = f-.1 


2 Tu 


4 Tu 


(37) 


,,  72+1  ,,  1 

4 = = 4+  — 

2 2 Tu  c 4TU 


(38) 


For  MSK-1,  the  modulation  functions  defined  for  the  time  interval  kTh  < t < (k+ 1)  T,  are 


3k(t)  = 4 


cos 


n t 


2 Tu 


(39) 


bit)  = - 0.sm 


n t 


2 T, 


rk^  = 1 


bit)  = - 


71  t 


2 T, 


hQk + Qi) 


(40) 

(41) 

(42) 


where  the  data  values,  dk,  are  used  to  create  two,  staggered,  x/2-rate  sequences  as  follows: 


1)  for  k even  and  0 < k < N 


h - = dk 


Qk  = dk~  i 

2)  for  k odd  and  1 < k < N 

Qk = Qk+i  = dk 


dk  - dk~  i 


Since  the  sequences  are  staggered,  only  one  of  the  sequences  can  change  its  value  at  each  data  epoch, 
kTb.  The  phase  functions  for  odd  and  even  k are, 


n t 


2 Tu 


dkdk- 1 + — [l“4-i]  for  k even 
2 


(43) 


(44) 


(45) 
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§ k(fi  = - dkdkl  + ^-[1-  dk]  for  k odd 
and  the  waveform  for  the  k-th  time  interval  is: 

lccos[27i/c*+(|)J 

^sin[2rc/e*+4)J 

= ikAcC0S  271 fy-  ^rJkQk+ 

= Accos2nfct--^IkQk+^[l-Ik]  + <\>c 


(46) 


(47) 


For  the  demodulation  of  MSK-1,  the  phase  differences  are  used  to  estimate  IkQk. 


A<t>*  = ~\hQk 

ikQk  = ~~^k 

71 

dk  = hQkdk-x  = --A<t>X-i 

71 


(48) 

(49) 

(50) 
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Figure  1 Basic  Quadrature  Detector 
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Figure  2 
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LPZF  I f Hilbert  Transform 


Figure  3 Generalized  Costas  Demodulator  (GCD) 
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